Microscopic analysis of the coherent optical generation and the decay of charge and 
spin currents in semiconductor heterostructures 
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The coherent optical injection and temporal decay of spin and charge currents in semiconductor 
heterostructures is described microscopically, including excitonic effects, carrier LO-phonon and 
carrier-carrier scattering, as well as nonperturbative light-field-induced intraband and interband 
excitations. A nonmonotonous dependence of the currents on the intensities of the laser beams is 
predicted. Enhanced damping of the spin current relative to the charge current is obtained as a 
consequence of spin-dependent Coulomb scattering. 

PACS numbers: 72.25.Fe, 42.65.-k, 72.25.Rb 
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One of the fundamental principles of quantum mechan- 
ics is that superpositions of wave functions lead to inter- 
ference phenomena which depend on the relative phase 
differences. In the coherent regime, e.g., shortly after an 
external perturbation has generated a nonequilibrium sit- 
uation, such quantum mechanical interference effects can 
be employed for the coherent control of dynamical pro- 
cesses in atomic^ molecular, biological, and semiconduc- 
tor systems 0, IE IE E IE IE] ■ Many of the measurements 
and proposals in this area make use of the coherent evo- 
lution of electronic excitations induced by specially de- 
signed optical laser pulses, e.g., sequences of phase-locked 
or suitably chirped beams. 

In semiconductors, the ultrafast coherent generation of 
photocurrents using two light beams with frequencies to 
and 2ui has attracted considerable attention 0, IE • As 
shown in Ref. flfl the same type of interference scheme 
can also be employed to create pure spin currents which 
are not accompanied by any charge current. Such spin 
currents generated on ultrafast time scales have been ob- 
served in semiconductors [ill and could be useful for 
future applications in the area of spintronics HEEl As 
shown recently, it is even possible to control photocur- 
rents via the carrier-envelope phase flE Il6l | which makes 
this scheme also interesting for optical metrology. Fur- 
thermore, for disordered semiconductors it has been pre- 
dicted that sequences of temporally delayed excitation 
pulses can be used to induce current echoes ^jj- 

The coherent generation of photocurrents in semicon- 
ductors by two light fields with frequencies to and 2u> 
satisfying 2Twj > £gap > Huj, where £gap is the band 
gap energy, has first been described in terms of nonlinear 
optical susceptibilities which have been obtained on the 
basis of band structure calculations [E ^3] • In this frame- 
work, the optically-induced intra- and interband transi- 
tions are treated within Fermi's golden rule. The inter- 
ference between intra- and interband excitations leads to 
electron and hole distributions which are not symmetric 
in k-space corresponding to a nonequilibrium situation 
with a finite current. The dynamics of the generation 



process has been analyzed using Bloch equations. This 
approach has been applied to disordered semiconductors 
within a two-band model |T^ | and to ordered quantum 
wells within a multiband formalism ^E] • Although the re- 
laxation of the photocurrent by carrier LO-phonon scat- 
tering bulk GaAs has been analyzed 20], in most of the 
existing publications the temporal decay of the charge 
and spin currents is still modeled by phenomenological 
decay times ^E EE EE E] ■ The Coulomb interaction 
among the photoexcited carriers has so far not been the 
focus of particular attention. However, excitonic effects 
in quantum wires have recently been addressed on the 
Hartree-Fock level |2lj . 

In this letter, we present and analyze a microscopic 
many-body theory that is capable of describing the dy- 
namical generation, the coherent evolution, and the de- 
cay of charge and spin currents. Our approach is based 
on the semiconductor Bloch equations (SBE), i.e., the 
equations of motion for the optical polarization and the 
carrier populations [22! l2.'"il| . The equations nonpertur- 
batively include the light-field-induced intraband and in- 
terband excitations without a rotating wave approxima- 
tion. The approach is thus capable of describing the non- 
trivial dependence of the initially generated currents on 
the intensities of the incident fields beyond the perturba- 
tive regime. The coherent part of the SBE contains the 
Coulomb interaction among the photoexcited carriers on 
the Hartree-Fock level, i.e., excitonic effects and Coulom- 
bic nonlinearities due to energy and field renormaliza- 
tions. As correlation contributions we include carrier- 
LO-phonon and carrier-carrier scattering at the second 
Born-Markov level EE El. 



The dynamical optoelectronic response is analyzed us- 
ing the Heisenberg equations of motion for the carrier 
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populations n% k = (aLk a CCT k) and «ak = 1 
and the interband polarization p CTk 
a Aerk ( a A<rk) creates (destroys) an electron with wave vec- 
tor k and spin a in band A. The resulting SBE including 
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where a = e, h, Vq denotes the Coulomb interaction po- 
tential, and d^. is the interband transition dipole. The 
terms given explicitely in Eqs. describe the dy- 

namics on the time-dependent Hartree-Fock level. The 
incoherent contributions are denoted by | coll . Here, these 
collision terms describe carrier LO-phonon and carrier- 
carrier scatterin g in the second-order Born-Markov ap- 
proximation |22ll23| . 

Except for the explicit inclusion of the spin, equations 
similar to Eqs. Q~© have been used to describe the 
optoelectronic response of semiconductor superlattices 
in the presence of static and Terahertz fields 0> E3- 
In these studies, the optical fields have been considered 
to generate interband transitions via E(t) ■ d™^ and the 
static and Terahertz fields lead to the intraband accel- 
eration via E(f) • V k . This distinction is useful if the 
involved fields are characterized by very different frequen- 
cies. However, since we here describe the generation of 
currents by the interference of optical fields with frequen- 
cies to and 2to this simplification is not possible and thus 
the total field is considered for both the intra- and inter- 
band excitations. 

Numerical solutions of Eqs. provide the time- 

dependent polarization and populations. The popula- 
tions in the valence and conduction bands determine the 
charge and spin current densities which are given by J = 
e E<xk v Xk " e E ff k v Xk and S = lE^k^Xk " 
I J2ak crv " 71 CTk' respectively, where v A = V k eA k /?i is the 
group velocity. 

The coherent generation of currents is due to ma- 
terial excitations which are not symmetric in k-space. 
In such situations, solutions of Eqs. CQ~© including 
the scattering contributions are quite demanding. In 
order to keep the numerical requirements within rea- 
sonable limits we consider for the present analysis two 
model systems: a one-dimensional (ID) quantum wire 
and a two-dimensional (2D) quantum well representa- 
tive of a GaAs/AlGaAs system. The electronic band 
structure is described in effective mass approximation us- 
ing m c =0.067 rrio an d m v = —0.457 too, and the band 
gap is -E gap =1.5 eV. The interband transition dipoles are 



taken as d| k = d cv (l,i,0) and dj k = d cv (l, — i,0) with 
d cv — 3 eA, i.e., we use the usual circularly polarized 
dipole matrix elements which describe heavy hole to con- 
duction band transitions in quantum wells close to the 
L-point [2^|. The incident laser field is given by 



E(t) = e " A v( e ~ 

1/—U),2U) 



C/^V^-^+ac.), (3) 



where e„ denotes the polarization, A v the amplitude, and 
4> v the phase of the field of frequency v = u> and 2w, 
respectively. Both frequency components are Gaussian 
shaped pulses with a duration determined by tl . For the 
case that both field components arc linearly polarized in 
x-direction, i.e., e w = e2 U = the photoexcitation pro- 
duces a pure charge current since the two spin systems 
are excited identically. For the case of linear perpendic- 
ularly polarized pulses, i.e., e w = and e 2u — e y , a 
pure spin current with no accompanying charge current 
is generated. Since the charge (spin) current is propor- 
tional to sia(<j> 2uj — 2^) (cos(0 2w — 20 w )) 01 we use 
<f> 2uJ — 2(f>u = tt/2 (<f>2ui ~ 20 w = 7r) in our calculations to 
obtain maximal currents. 

Time-dependent electron and hole distributions of a 
quantum well in k-space are shown in Fig. ^ The car- 
riers are initially generated with a combined excess en- 
ergy of 150 meV above the band gap. Due their smaller 
mass about 130 meV of the kinetic energy is given to 
the electrons and the excess energy of the holes is only 
about 20 meV. Therefore, the electron relaxation is con- 
siderably influenced by LO-phonon emission whereas the 
hole distribution is only weakly affected. Immediately af- 
ter the excitation, the electron and hole distributions are 
very similar, see Fig.^ Both are generated nonuniformly 
on a ring with radius f=a 5/ao. Due to the quantum inter- 
ference, the distributions are larger for positive k x than 
for negative k x . Thus the distributions have a nonvan- 
ishing positive average momentum which corresponds to 
a current in x-direction. In the course of time, the dis- 
tributions relax towards quasi-equilibrium distributions. 
For long times, due to their larger mass the distribution 
of the holes is wider than that of the electrons. 

Figure Ufa) demonstrates that for the considered ex- 
citation conditions both carrier LO-phonon and carrier- 
carrier scattering contribute significantly to the current 
dynamics. If only carrier LO-phonon scattering is taken 
into account, the charge and spin currents decay simi- 
larly. This decay is not exponential, however, its onset 
can be described by an exponential decay with time con- 
stant 240 fs. When carrier-carrier scattering is included, 
the currents decay more rapidly and, in particular, the 
spin current decays faster than the charge current. Addi- 
tional calculations which omit specific contributions have 
revealed that this difference is due to Coulomb scatter- 
ing between carriers with different spin. When a charge 
current is excited (e w = e-i^ — e x ), the electron distri- 
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FIG. 1: (color) Left (right) column: Contour plots of the 
electron (hole) distributions of a quantum well in k-space at 
t =50, 100, 150, and 400 fs (from top to bottom), respec- 
tively. The incident pulses have a duration of tx=20 fs, the 
amplitudes are A u = 2^2^ = 108Ao, with Ao = Eo/eao ~ 
4 x 10 3 V/cm, where Eq is the three-dimensional exciton Ry- 
dberg, and 2hu> = 1.65 eV. The density of the photoinjected 
carriers is N = 10 11 cm" 2 and the temperature is T — 50 K. 
ao is the three-dimensional exciton Bohr radius. 



in this case both spins have the same nonvanishing av- 
erage momentum. Since carrier-carrer scattering only 
exchanges momentum among the carriers, this average 
momentum is not reduced by the Coulomb scattering. 
The situation is, however, different when a spin current 
is excited (e w = e x , = e y)- In this case the electron 
distributions satisfy 7i| k = n% where k has the same y- 
component as k but its negative x-component, i.e., the 
average momenta of the two spins point into opposite 
directions. Therefore, the total electron momentum van- 
ishes and Coulomb scattering leads to a decay of the av- 
erage momenta of the spin-up and spin-down electrons. 
Fig. Hfb) shows that qualitatively similar results are ob- 
tained for quantum wires. However, due to the smaller 
ID phase space, the scattering is reduced and thus the 
decay times are longer than in 2D. The results presented 
in Fig. 13 and of additional evaluations show that i) at 
low densities the charge and spin currents decay with the 
same time constant due to carrier LO-phonon scattering 
and ii) with increasing density, carrier-carrier scattering 
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FIG. 2: (a) Time-dependent charge (solid) and spin (dash- 
dot) currents of a quantum well for the same parameters as 
in Fig. Also shown is the identical decay of both currents 
if only carrier LO-phonon scattering is considered (dashed). 
The thin solid lines represent exponential decays with time 
constants of 240, 155, and 125 fs, respectively, (b) Same as 
(a) for a quantum wire. The density of the photoinjected 
carriers is N = 5 x 10 5 cm" 1 and the other parameters are 
the same as in (a) . The thin solid lines represent exponential 
decays with time constants of r = 1250, 900, and 740 fs, 
respectively. 



speeds up the decay of both currents, however, the rate of 
change is larger for the spin than for the charge current. 

Since 2D calculations are numerically extensive, we 
present in the following results obtained in ID. Qual- 
itatively similar effects should also be present in 2D. 
Fig. |3f a) shows how the injected carrier density depends 
on the field amplitude A — A 2 ^ — A u /2. In the con- 
sidered regime, the density increases approximately lin- 
ear with A. However, the current depends on A in a 
strongly nonlinear fashion, see Fig.[3{b). For small am- 
plitudes, the nonsymmetric k-space carrier distributions 
increase with A without significant distortion. There- 
fore, in this regime both the total density and the charge 
current become larger with increasing A. At a certain 
excitation level the peak of the generated distribution be- 
comes comparable to 1 and further excitation is reduced 
by the Pauli principle. Thus, in the high-field regime the 
asymmetry of the k-space distribution and the current 
decrease with increasing field amplitude. 

Figures |3{c) and (d) show how the carrier and the 
charge current densities depend on the ratio of the field 
amplitudes x = A^jAi^ for a fixed amplitude of the 
2uj field A 2lu = 128A a . At x = 0, only the 2uj field is 
present which generates via interband excitations sym- 
metric carrier distributions in k-space, i.e., no current. 
A finite current is present only if x is finite since the 
combined action of inter- and intraband excitations is 
required for the generation of nonsymmetric k-space car- 
rier distributions. In the limit of small x the current 
increases as x 2 in agreement with a perturbative analysis 
of the light-matter coupling. In this regime, the carrier 
density increases only slightly with x since a weak w field 
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FIG. 3: Dependence of (a) the carrier density and (b) the 
charge current density at different times on the amplitude 
A = A2u> = Au,/2 of the incident pulses in a quantum 
wire. Dependence of (c) the carrier density and (d) the 
charge current density on the amplitude ratio x = A^jA^ 
for A^ = 128 A . Here, we use 2hoj = 1.62 eV and T = 50 K. 



predominantly redistributes the carriers in k-space by in- 
troducing asymmetries. Beyond this perturbative regime 
we find an interesting dependence of both the carrier den- 
sity and the current on x. The largest current is obtained 
for x « 4, i.e., k, 512Ao. It can be expected that 
for optical frequencies, the x which generates the largest 
current is always larger than 1, since intraband excita- 
tions are relatively smaller than interband excitations. 
We find that the optimal value of x increases with de- 
creasing . Thus in the limit of weak fields A u should 
be chosen much larger than A2 U if one is interested in 
generating large currents. Our calculations have shown 
that the intensity dependence of the spin current is very 
similar that of the charge current, see Fig- EI 

In summary, the coherent optical injection and the de- 
cay of charge and spin currents in semiconductor het- 
erostructures is described by a microscopic many-body 
theory. We find that due to Coulomb scattering between 
carriers of different spin at elevated excitation levels the 
spin current decays more rapidly than the charge cur- 



rent. An interesting nonmonotonous dependence of the 
currents on the intensities of the laser beams is predicted. 
These results should stimulate further experimental re- 
search in this direction. 

This work is supported by the Deutsche Forschungs- 
gemeinschaft (DFG), by the Ministry of Education and 
Research (BMBF), and by the Center for Optodynam- 
ics, Philipps University, Marburg, Germany. We thank 
the John von Neumann Institut fur Computing (NIC), 
Forschungszentrum Julich, Germany, for grants for ex- 
tended CPU time on their supercomputer systems. 



[1 
[2 

[3. 

K 

[5. 
[6 

[-: 

[8 
[9 
[10 

[11 
[12 

[13 

[14 

[is; 

[16 
[17 

[is; 

[19 

[20 
[21 

[22 

[23 

[24 
[25 
[26 



P. Brumer and M. Shapiro, Chem. Phys. Lett. 126, 541 
(1986). 

W.S. Warren, H. Rabitz, and M. Dahlem, Science 259, 
1581 (1993). 

R.J. Gordon and S.A. Rice, Ann. Rev. Phys. Chem. 48, 
601 (1997). 

T. Brabec and F. Krausz, Rev. Mod. Phys. 72, 545 
(2000). 

J.L. Herek et al, Nature 417, 533 (2002). 

M. Shapiro and P. Brumer, Rep. Prog. Phys. 66, 859 

(2003). 

E. Dupont et al, Phys. Rev. Lett. 74, 3596 (1995). 
R. Atanasov et al, Phys. Rev. Lett. 76, 1703 (1996). 
A. Hache et al, Phys. Rev. Lett. 78, 306 (1997). 
R.D.R. Bhat and J.E. Sipe, Phys. Rev. Lett. 85, 5432 
(2000). 

M.J. Stevens et al, Phys. Rev. Lett. 90, 136603 (2003). 
J. Hiibner et al, Phys. Rev. Lett. 90, 216601 (2003). 
S.A. Wolf et al, Science 294, 1488 (2001). 
I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

T.M. Fortier et al, Phys. Rev. Lett. 92, 147403 (2004). 
PA. Ross et al, J. Opt. Soc. Am. B 22, 362 (2005). 

C. Schlichenmaier et al, Phys. Rev. B 65, 085306 (2002). 
A. Najmaie, R.D.R. Bhat, and J.E. Sipe, Phys. Rev. B 
68, 165348 (2003). 

D. H. Marti, M.-A. Dupertuis, and B. Deveaud, Phys. 
Rev. B 69, 035335 (2004). 

P. Krai and J.E. Sipe, Phys. Rev. B 61, 5381 (2000). 
D.H. Marti, Ph.D. thesis, Ecole Polytechnique Federale 
de Lausanne, Switzerland, 2003. 

H. Haug and S. W. Koch, Quantum Theory of the Opti- 
cal and Electronic Properties of Semiconductors, 4th ed. 
(World Scientific, Singapore, 2004). 

W. Schafer and M. Wegener, Semiconductor Optics and 
Transport Phenomena, (Springer, Berlin, 2002). 
T. Meier et al, Phys. Rev. Lett. 73, 902, (1994). 
T. Meier et al, Phys. Rev. Lett. 75, 2558 (1995). 
These matrix elements are clearly not realistic for quan- 
tum wires. This model, however, allows us to perform de- 
tailed investigations of the currents in a ID system with 
rather modest numerical requirements. The obtained re- 
sults are useful since the fundamental dependencies and 
the differences between the decay of the charge and spin 
currents are similar in 2D. 



